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Heavy fermion (HF) materials exhibit a rich array of phenomena due to the strong Kondo cou¬ 
pling between their localized moments and itinerant electrons. A central question in their study is 
to understand the interplay between magnetic order and charge transport, and its role in stabiliz¬ 
ing new quantum phases of matter. Particularly promising in this regard is a family of tetragonal 
intermetallic compounds CeTX 2 {T = transition metal, X — pnictogen), that includes a variety of 
HF compounds showing T-linear electronic specific heat Ce ~ 7T, with 7 ~ 20-500 mJ-mol~^ K“^, 
reflecting an effective mass enhancement ranging from small to modest. Here, we study the low- 
temperature field-tuned phase diagram of high-quality CeAgBi 2 using magnetometry and trans¬ 
port measurements. We find an antiferromagnetic transition at Tn = 6.4 K with weak magnetic 
anisotropy and the easy axis along the c-axis, similar to previous reports (Tjv = 6.1 K). This scenario, 
along with the presence of two anisotropic Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions, 
leads to a rich field-tuned magnetic phase diagram, consisting of five metamagnetic transitions of 
both first and second order. In addition, we unveil an anomalous Hall contribution for fields H < 54 
kOe which is drastically altered when H is tuned through a trio of transitions at 57, 78, and 84 kOe, 
suggesting that the Fermi surface is reconstructed in a subset of the metamagnetic transitions. 


In heavy fermion (HF) materials, the Kondo coupling 
between local moments and itinerant electrons plays a 
central role in determining magnetic and transport prop¬ 
erties, particularly at low temperatures. Classic exam¬ 
ples of the unusual behavior include quantum critical¬ 
ity in YbRh 2 Si 2 [1], unconventional superconductivity in 
CeCoIns [2], and metamagnetism in CeRu 2 Si 2 [3]. Ce- 
based HF materials which host such exciting properties 
often crystallize in tetragonal structures and their ground 
state is on the border of antiferromagnetism (AFM). 

Our focus in this paper is the HF family CeTX 2 {T = 
transition metal, X = pnictogen), a class of intermetallic 
compounds that crystallize in the tetragonal ZrCuSi 2 - 
type structure (space group PA/nmm) with a stacking 
arrangement of CeX-T-CeX-X layers. Great attention 
has been given to the antimonide {X = Sb) members due 
to the presence of anomalous ferromagnetism in CeAgSb 2 
and, more recently, due to the report of field-induced 
quantum criticality in CeAuSb 2 [4, 5]. Although the in¬ 
vestigation of the antimonides is abundant, fewer reports 
can be found on the bismuthide (A = Hi) members. Re¬ 
cent studies of high-quality CeCuBi 2 and CeAuBi 2 re¬ 
vealed AFM ordering temperatures of Tjv =16 K and 
12 K, respectively. Both materials show a moderate mass 
enhancement from T-linear electronic specific heat (7 on 
the order of 100 mJ-mol“^ K“^) and exhibit a single spin- 
flop transition with field applied along the c-axis [6, 7]. 

CeAgBi 2 , a third member of the isovalent series, has 
been reported previously to order antiferromagnetically 
(Tjv = 6.1 K) at zero magnetic field and to undergo 
three field-induced magnetic transitions at 2 K (iJ^ = 
35, 50,83 kOe) [8, 9]. However, a detailed analysis of the 
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underlying interactions and a deep understanding of the 
phase diagram at low temperatures is still missing. 

Single crystals of CeAgBi 2 were grown from Bi-flux 
with starting composition Ce:Ag:Bi=l:a;„omma;:8 (1 < 
^nominal < 3). The mixture was placed in an alumina 
crucible and sealed in a quartz tube under vacuum. The 
sealed tube was heated up to 1050°C for 8 h and then 
cooled down at 10°C/h. The excess of Bi flux was re¬ 
moved by centrifugation after 24 h of annealing at 500° C. 
Single crystals with dimensions ~ 3 x 3 x 0.5 mm^ were 
ground and their crystal structure was checked by X-ray 
powder diffraction experiments using Cu Ka radiation at 
room temperature. Several single crystals from different 
batches were also submitted to elemental analysis using 
a commercial Energy Dispersive Spectroscopy (EDS) mi¬ 
croprobe coupled to a FEG SEM microscope. From the 
EDS analysis, we have extracted the actual XAg concen¬ 
tration. The precision of the analysis was calculated by 
a/^/N, where a is the standard deviation of the mea¬ 
surements, and N is the number of points analyzed. 

Magnetization measurements were performed using a 
commercial superconducting quantum interference de¬ 
vice (SQUID) down to 1.8 K. Below 1.8 K, cantilever- 
based torque magnetometry was used to measure the 
magnetization in a dilution refrigerator. Electrical re¬ 
sistivity measurements were made using a low-frequency 
ac resistance bridge and a four-point configuration. 
Hall data was obtained by measuring in both positive 
and negative applied fields and taking the difference 
{[R{H-\-) — R{H—)] /2). Due to hysteresis in the tran¬ 
sitions data was paired for subtraction based on whether 
the magnitude of the field was increasing or decreasing. 
All measurements reported here, except for some of the 
magnetic susceptibility measurements, were made with 
the applied field parallel to the c-axis. For resistivity 
measurements, the current was applied in the ab-plane. 



2 


I. EXPERIMENTAL RESULTS AND CEF FITS 

Fig. la shows the temperature dependence of in-plane 
resistivity, p{T), down to 0.5 K at zero magnetic field. 
At high temperatures {T > 200 K), p{T) shows metallic 


a 



FIG. 1: (a) Temperature dependence of in-plane resistivity, 

p(T), in zero applied field. Top inset shows the low-T region 
where a kink is observed at Tm- Bottom inset shows Tjv as 
a function of Ag occupancy measured by EDS. (b) Inverse 
susceptibility data for fields parallel (x||) and perpendicular 
(X-l) to the c-axis. Inset shows low temperature x(T) data. 
The solid lines are fits to the data using the model described 
in the text, (c) Magnetic contribution to specific heat in zero 
field. Top inset shows the integration of specific heat. Snorm 
is defined as S/(J?ln2). Bottom inset shows p{T) in three 
different magnetic fields (in kOe). 


behavior, decreasing linearly with decreasing tempera¬ 
ture. However, further decrease in temperature reveals a 
resistivity minimum, followed by a logarithmic increase 
due to incoherent Kondo scattering. Below ^ 25 K, p{T) 
drops abruptly, suggesting that this is the energy scale of 
either CEF depopulation or Kondo coherence. We will 
discuss these possibilities below. The kink in resistiv¬ 
ity observed at T^r = 6.4 K, shown in the top inset of 
Fig. la, indicates the transition to the AFM phase. This 
represents a 0.3 K higher ordering temperature compared 
to previous studies of CeAgBi 2 [8, 9]. Although this 
increase in Tjv is somewhat small, EDS measurements 
reveal that it is caused by a substantial decrease both 
in the number of vacancies and in the inhomogeneity at 
the Ag site. The bottom inset of Fig. la shows the lin¬ 
ear dependence of T/v on the occupation at the Ag site, 
a^EDS- For the most deficient samples, with Ag occupa¬ 
tion of 82(4)%, the transition temperature matches previ¬ 
ous studies, = 6.1 K. For the best samples obtained 
to date, the Ag occupancy reaches 87(2)%, confirming 
the trend that transition-metal deficiency is an intrinsic 
feature of this family of compounds. We note that less 
deficient samples are also accompanied by higher resis¬ 
tance ratios and lower residual resistivity. 

Like other Ce-based bismuthides, CeAgBi 2 also ex¬ 
hibits magnetic anisotropy. Fig. lb shows the temper¬ 
ature dependence of the inverse magnetic susceptibility, 
l/x(r), when a magnetic field of = 1 kOe is applied 
parallel (yy) and perpendicular (x_l) to the crystallo¬ 
graphic c-axis. The inset of Fig. lb presents the low- 
temperature x(T) data in which a sharp peak is observed 
at Tpj. The ratio X||/X-L ~ 3-5 at T^r is mainly deter¬ 
mined by the tetragonal CEF splitting and reflects the 
low-T Ce^“'" single ion anisotropy. This ratio is smaller 
than what is found in other bismuthides, suggesting a 
smaller CEF splitting between the ground state and the 
first excited state as well as a less anisotropic ground 
state. This scenario will be confirmed below. 

At high temperatures (T > 150 K), x(F) is well- 
described by a Curie-Weiss (CW) law plus a T- 
independent Pauli term, x(F) = Xo + C/{T — Ocw)- 
We obtain an effective moment of /ie// ~ 2.5(l)/iB for 
both directions and also for the polycrystalline average, 
in agreement with the theoretical value of Peff ~ 2.54/iB 
for Ce^’*' free ions. On the other hand, the Ocw values 
are anisotropic, with = 5.7 K and 0± = —8.5 K. For 
the polycrystalline averaged data (not shown), we obtain 
Op = —4 K, consistent with AFM order at ~ 6 K. 

We now further explore the role of anisotropic inter¬ 
actions and CEF effects in the magnetic properties of 
CeAgBi 2 . To this end, we analyze the experimental data 
using a mean field model including two mean anisotropic 
interactions (zafM'^afm and zfM’^fm) between nearest- 
neighbors, which do not contain directional information. 
We also take into account the tetragonal CEF Hamil¬ 
tonian Hcef = B 2 O 2 + + BjOj, where z is the 

number of nearest neighbors, are the CEF param¬ 
eters, and O" are the Stevens equivalent operators ob- 
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CEF parameters 


^2“ 

-1.78 K 


B'i 

0.168 K 


Bt 

0.71e-3 K 


zfmJfm 
-0.89 K 


zafmJafm 

1.35 K 


Energy levels and wave functions 


E{K) 

1 - 5/2) 

1 - 3/2) 

1 - 1/2) 

|±l/2) 

1 ± 3/2) 

1 ± 5/2) 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

19 

1 

0 

0 

0 

0 

0 

19 

0 

0 

0 

0 

0 

-1 

61 

0 

0 

-1 

0 

0 

0 

61 

0 

0 

0 

-1 

0 

0 


TABLE I: CEF parameters, energy levels, and wave fnnctions of CeAgBi 2 single crystals obtained from the best fits of the 


magnetic susceptibility data to the model described in the text. 


tained from the angular momentum operators. A more 
detailed description of the model can be found in Ref. 
[ 10 ]. 

This model was used to simultaneously fit x(T), M{H) 
and Cmag{T)/T data in the entire range of tempera¬ 
ture. The best fit which reproduces the anisotropic 
susceptibility is shown by solid lines in Fig. lb and 
the extracted parameters for the CEF scheme and the 
exchange interactions are given in Table L We find a 
|J = 5/2, Jz = ±3/2) Kramers doublet ground state, 
separated from other excited doublets by 19 K and 
60 K. We also find, in fact, (i) the CEF splitting is 
much smaller as compared to other bismuthides and (ii) 
the ground state is mainly |J = 5/2, = ±3/2) in¬ 
stead of the dominantly |J = 5/2, = ±5/2) found 

in CeCuBi 2 and CeAuBi 2 . Further, the dominant CEF 
parameter i ?2 obtained from the fits is similar to the 
value obtained using the high-temperature expansion of 
X: Bl = 10(6»j.-6»||)/[3(2J-l)(2J-b3)] =-1.48 K.[11] 
This result suggests that the effects of anisotropic inter¬ 
actions at high temperatures are smaller in CeAgBi 2 than 
in CeCuBi 2 and CeAuBi 2 . We also note that although 
the magnetic anisotropy and Tjq along the c-axis are well 
reproduced by our simple model, the single ion CEF ef¬ 
fect is not able to capture all held-induced transition in 
M{H) data at 1.8 K. 

Fig. Ic shows the temperature dependence of the spe¬ 
cific heat of CeAgBi 2 in zero field. To determine the 
magnetic contribution to specific heat, the data obtained 
from non-magnetic reference compound LaAgBi 2 was 
subtracted from the result. There is a single peak cor¬ 
responding to Tat. We note that this sample is from 
an earlier batch, so it has higher silver vacancies and a 
lower Tjv. The top inset of Fig. Ic shows the integra¬ 
tion of specific heat over temperature, i.e., the recovered 
magnetic entropy. At T^r, the recovered entropy is only 
about 80% the ground-state doublet (Rln2). We note, 
however, that there is a broad feature centered around 
8 K, which is consistent with a Schottky anomaly gen¬ 
erated by the CEF splitting of 19 K. This result indi¬ 
cates that the first excited CEF state is already par¬ 
tially occupied at T^r, giving an additional entropy to 


the expected Rln2. Therefore, the reduction of entropy 
is larger than 20%, although a precise calculation cannot 
be evaluated. It is unlikely that a reduction of entropy 
of more than 20% can be fully accounted for by classical 
magnetic fluctuations. For example, numerical calcula¬ 
tions of the Ising model on a simple cubic lattice find an 
entropy reduction of just under 20% at Tat, with other 
three dimensional geometries showing an even smaller 
reduction. [12] We thus attribute this reduction to both 
magnetic frustration and hybridization between Ce 4/- 
electrons and conduction electrons. This conclusion is 
motivated by the enhanced effect masses found in dHvA 
experiments [8] and the magnetic frustration that often 
arises from competing exchange interactions. 

Due to the similar energy scales of the first excited 
CEF state (19 K) and the peak in resistance (24 K), we 
performed high-temperature resistance measurements in 
several different helds. The results are presented in inset 
of Fig. Ic. If the peak near 24 K were due to crystal 
field depopulation, then increasing the field should lead 
to an increase of the energy difference between the low¬ 
est Zeeman split state of the ground state and the upper 
Zeeman split state of the first excited CEF state. This in 
turn would lead to an increase of the temperature of the 
resistance peak as the field is increased, which is in fact 
observed experimentally. We note that, in heavy fermion 
compounds where the energy difference between Kondo 
coherence and Acef is small, the onset of Kondo coher¬ 
ence (T*) usually occurs at temperatures of the order 
of Acef/ 2. This suggests that T* ~ 10 K in CeAgBi 2 . 
However, the experimental lack of two distinct resistance 
maxima as is typically observed when the energy differ¬ 
ence is large[13] leads to an ambiguity as to whether the 
resistance peak is due to CEF splitting of Kondo coher¬ 
ence. 

We now turn to the analysis of the low-temperature 
phase diagram as a function of an applied magnetic field 
along the c-axis. Fig. 2a shows magnetization and (trans¬ 
verse) magnetoresistance (MR) of CeAgBi 2 at 100 mK. 
As the field is swept from 0 to 90 kOe, we find a se¬ 
quence of field-induced metamagnetic transitions near 
H = 34,37,54, 78, and 84 kOe, indicating a complex 
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phase diagram for T < T^r. Above the 84 kOe transi¬ 
tion, the magnetization saturates just below 1.9 fis/Ce, 
which is slightly lower than the saturation value of 2.1 
/rs/Ce previously reported, and may be due to a non¬ 
linearity between the deflection of the cantilever and the 
corresponding change in capacitance. Such large satu¬ 
rated value is somewhat surprising since the largest pos¬ 
sible c-axis magnetization of a | ±3/2) CEF ground state 
is (3/2)5i./rB, i.e., 1.29 /is/Ce. Hence, our results indi¬ 
cate that the Zeeman effect induces a change in the CEF 
ground state from | ± 3/2) to | ± 5/2) when H ^ 84 kOe. 




FIG. 2: Field-tuned magnetic and transport properties of 
CeAgBi 2 . (a) (blue) and magnetization (green) versus 
field at 100 mK. p^^x data is for both increasing and decreas¬ 
ing fields, showing hysteresis for transitions centered near 54 
and 78 kOe. (b) Magnetization versus field for decreasing 
field-sweep at temperatures 500 mK and lower. Sweep rate 
was 500 Oe/min. Inset shows hysteresis in magnetization de¬ 
pending on field sweep direction. 

Of the observed transitions, the one near 78 kOe ap¬ 
pears to be strongly first-order as shown by magnetiza¬ 
tion in Fig. 2b. As T is lowered from 500 mK to 35 mK 
a sharpening of the transition is observed and the tran¬ 
sition eventually becomes step-like. Field sweeps per¬ 
formed around this transition show signs of irreversibility, 
namely a ~ 4 kOe wide hysteresis loop, providing further 
evidence of first-order behavior. This in contrast with 
other Ce-based bismuthides. For instance, other mem¬ 


bers of the isovalent series such as CeCuBi 2 (T^y = 16 K) 
and CeAuBi 2 (T^y = 12 K) show only a single clear spin- 
flop transition. We note, however, that multiple steps 
have been observed in previous reports on more deficient 
CeCuBi 2 {Tjsr = 11 K) and even CeCuBi 2 (T^r = 16 K) 
shows hints to a second transition just before the main 
one. It is also noteworthy that hysteresis loops show a 
significant difference (~ 2 kOe) at the transition centered 
near 54 kOe as well. The change in magnetization, how¬ 
ever, is not nearly as abrupt as found in the 78 kOe tran¬ 
sition. It is well-known that transport measurements in 
real materials are complicated by disorder. Thus, as this 
transition is clearly not strongly first-order, the slight 
amount of irreversibility is likely attributed to pinning 
due to crystallographic defects (e.g., inherent silver de¬ 
ficiency). In fact, a small amount of hysteresis (on the 
order of a hundred gauss) was also observed in magneti¬ 
zation and transport at the remaining three transitions. 

The rich phase structure in CeAgBi 2 is likely due to 
the weak anisotropy combined with anisotropic exchange 
parameters with opposite signs, which lead to magnetic 
frustration. Further, the close energy scales of the CEF 
splitting (^ 19 K), Kondo coherence (~ 10 K) and the 
AFM {Tn ~ 6.4 K) generates a complex response of the 
physical properties to the application of magnetic field. 

The MR also reflects the multi-step phase structure, 
tracking each metamagnetic transition with a sharp step, 
indicating the presence Kondo coupling between the Ce 
magnetic moments and itinerant p-electrons from Bi. 
The overall trend of MR is linear, a feature shared with 
the (non-magnetic) compounds LaAgBi 2 and LaAgSb 2 , 
where it is attributed to an underlying Dirac dispersion 
for Bi/Sb itinerant electrons [14, 15]. A preliminary ab 
initio calculation [16] for CeAgBi 2 indicates that the Bi 
p-electrons also have a Dirac dispersion, suggesting a sim¬ 
ilar origin for the underlying linear MR. However, we note 
an alternative explanation for the increase in MR is an 
enhancement in spin-disorder scattering in AFM materi¬ 
als in an applied field [17, 18]. Note that the MR exhibits 
a sharp drop and hysteretic behavior near the 79 kOe 
transition, lending further support that it is first-order. 

Motivated by this evidence for coupling of itinerant 
electrons and local moments, we performed measure¬ 
ments of the Hall resistivity pxy to further elucidate the 
nature of the different phases. Like the MR, the Hall 
effect tracks each metamagnetic transition (Fig. 3). 

The decrease in p^y roughly parallels the increase in 
magnetization over the first two transitions, but then 
sharply deviates from this trend as the field is increased 
across the trio of transitions at 57, 78, and 84 kOe, even 
changing its sign between ~ 60 — 80 kOe. Further insight 
is afforded by fitting the Hall resistivity to the standard 
form pxy{H) = RhH RmM{H) that includes Hall ef¬ 
fect contributions from both the applied field and the in¬ 
duced magnetization. Using the data up to < 38 kOe 
to obtain the fit parameters, we find that this model fits 
the data extremely well for low fields H < 54 kOe, where¬ 
upon the measured pxy strongly deviates from the ex- 
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FIG. 3: Hall resistance {pxy) versus field, (a) Comparison of 
sweeping the field up vs. down at 150 mK. The dotted line is 
a least-squares fit to Pxy{H) — RhH + RmM{H). The fit was 
made to the region of zero field up to 37 kOe (green arrow) 
and extrapolated beyond that point, (b) In a second sample, 
comparison at different T while sweeping the field downward. 


pectation based on the model. While there is still some 
discrepancy between the fit and the data after the last 
transition {H > 84 kOe), this is relatively small com¬ 
pared to the much larger deviation between fit and data 
in the regime 54 kOe H < 84: kOe. As the tempera¬ 
ture is raised, the magnitude of the anomaly in p^y near 
54 kOe is decreased, as shown in Fig. 3b. At present, 
we do not have a complete understanding of the origin 
of this Hall effect anomaly, but we comment on some 
possible scenarios below. We note that while there is an 
additional deviation from linear behavior of pxy above 
^110 kOe, a previous study on lower-quality samples 
using pulsed fields up to 400 kOe [8] appears to rule out 
an additional transition. 


II. DISCUSSION 

Our results can be summarized in the form of a phase 
diagram of CeAgBi 2 with temperature and applied field, 
as shown in Fig. 4. The second T = 0 phase (for 
34 kOe < H < 37 kOe) is not seen above the temper¬ 
ature T > 2.5 K, illustrating the fragility of this phase 
to thermal fluctuations. Note the strong agreement be¬ 
tween the probes of specific heat, transport, and mag¬ 
netic structure when the corresponding measurements 
overlap. Of the sequence of transitions between the differ¬ 
ent low-temperature phases, the penultimate transition 
near 79 kOe appears to be strongly first-order, whereas 
the remaining transitions appear continuous. As noted, 
the Hall effect is well-described by a simple model includ¬ 
ing effects of both the applied field and induced magneti¬ 
zation very well below 54 kOe and reasonably well above 
84 kOe, but exhibits a strong deviation from this simple 
behavior in the intermediate field regime. 



FIG. 4: (a) Low temperature H-T phase diagram from 

GeAgBi2. Phase boundary was determined by resistivity, 
magnetization, and heat capacity measurements. 

Heavy fermion compounds CeAuSb 2 and YbAgGe 
both have some commonalities in the low-temperature 
phase diagram. For instance, CeAuSb 2 orders AF and 
has at least one magnetic transition before the suppres¬ 
sion of AF order to 0 K at 54 kOe.[5] Recent work sug¬ 
gests that the first transition might bifurcate into two 
at lower temperatures, similar to CeAgBi2.[19] YbAgGe 
also has a similarly complicated phase diagram as a func¬ 
tion of field which contains a region partially bordered by 
a first order transition [20]. Interestingly, in this region 
YbAgGe exhibits anomalous behavior attributed to the 
influence of a quantum bicritical point [21]. Unlike in 
YbAgGe, however, the set of transitions is only observed 
in CeAgBi 2 for H\\c. Further, the large anomalies in the 
Hall resistivity appear to be a feature unique to CeAgBi 2 . 

We now comment on potential explanations for the 
physics at hand. First, we attribute the sequence of 
metamagnetic transitions to the frustration from com- 
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peting anisotropic exchange interactions, as evinced by 
the susceptibility data. As in many HFs, the Hall data is 
more challenging to describe [22]. Its evident sensitivity 
to the magnetization suggests a significant contribution 
from the induced magnetization at low fields, yet a simple 
phenomenological model incorporating this contribution 
seems unable to describe the behavior at intermediate 
fields. Here we note that quantum oscillation measure¬ 
ments reveal electron effective masses in the range of 5 
to 7 me, suggesting heavy fermion physics plays a role in 
transport [8]; this could perhaps explain the anomaly, e.g. 
via a Fermi surface reconstruction. To explore this sce¬ 
nario, one would have to investigate, for instance, band- 
structure calculations in the presence of a magnetic field. 
We anticipate that future theoretical analysis, incorpo¬ 
rating details of the conduction electron band structure, 
as well as a careful treatment of their coupling to local 
moments, will shed more light on these issues [16]. 

An alternative scenario involves the ‘geometrical Hall 
effect’ from a non-collinear spin texture; while more 
exotic, such physics is now quite well-established as 
a consequence of Dzyaloshinskii-Moriya interactions in 
ferromagnets and has been also proposed to occur in 
anisotropic exchange AFMs, and recently also in heavy- 
fermion metals [23]. In this case, an anomalous Hall 
effect could be observed as a result of a non-vanishing 
Berry-phase curvature [24]. To this end, neutron diffrac¬ 
tion/scattering experiments would be valuable to unveil, 
for instance, the field regions where non-collinear spin 
structures are realized. 


III. SPIN MODEL AND METAMAGNETIC 
TRANSITIONS 


As mentioned above, our previous model was unable 
to capture all field-induced magnetic transitions. Finally, 
in the present section we introduce a spin model for the 
Ce local moments to discuss their field-induced response. 
Although more experimental investigation is clearly es¬ 
sential in order to clarify what type of AFM order is sta¬ 
bilized in CeAgBi 2 in each field region, here we sketch one 
possible scenario built on existing experimental data on 
the closely-related compound CeCuBi 2 , where the zero- 
field magnetic structure has already been well-established 
experimentally [7]. We defer a detailed analysis of the 
spin model to future work, and here only discuss a sim¬ 
ple mean-field picture of the magnetic order. 

For CeCuBi 2 , X-ray magnetic diffraction reveals that 
the Ce local moments order antiferromagnetically with 
an ordering wave vector (0,0,0.5) [7]. This corre¬ 
sponds to “up-up-down-down” magnetic ordering along 
the crystallographic c-axis and ferromagnetic (FM) or¬ 
dering across the basal a-b plane. (See Fig. 5). The mag¬ 
netic susceptibility data shows a large anisotropy consis¬ 
tent with a strong easy-axis along the c-axis. The appli¬ 
cation of a c-axis magnetic field results in a single spin- 
flop transition at a field of ~ 55kOe; for fields in the ab 


plane, no such transition was observed. Within a clas¬ 
sical approximation, a minimal spin model that may be 
inferred from these experimental results incorporates a 
single-ion anisotropy term and dominant spin exchanges 
between neighboring sites: 

h = J2 -s.-Y^h-s, -AYis^^f■ ( 1 ) 

ij i i 


Here, h is the (applied) magnetic field and A > 0 is the 
magnitude of an easy-axis single-ion anisotropy along the 
c-axis. Jij represents competing spin exchanges between 
neighboring sites; FM and AFM (or AFM and FM) spin 
exchanges for two symmetry distinct neighbors along c 
axis are denoted J and J', respectively, and a FM spin 
interaction in the a-b plane is denoted as J_l, as seen in 
Fig. 5. The competition between alternating ferro- and 
antiferromagnetic spin exchanges along the c-axis stabi¬ 
lizes the “up-up-down-down” magnetic structure along 
c-axis, with a spin polarized state FM on each a-b plane. 
In the limit A 3> Jij-, the moments behave like c-axis 
Ising spins, and the presence of a c-axis field leads to a 
direct magnetic transition from AF ordering to a fully- 
spin-polarized, paramagnetic state. However, a smaller 
single ion anisotropy A < Jij admits a spin-flop tran¬ 
sition at finite c-axis field, resulting in AF spin canting. 
This is consistent with the experimental data. Given that 
the minimal picture of competing spin exchanges and sin¬ 
gle ion anisotropy provide a reasonable explanation of 
magnetic ordering and field effect of CeCuBi 2 material, 
we now turn to a similar analysis of CeAgBi 2 . 




a 



FIG. 5: Lattice structure of CeAgBi 2 : Red, blue and green 
spheres represent Ce, Ag and Bi ions respectively and solid 
line shows one unit-cell, “up-up-down-down” like ordering of 
Ce moments is also shown along c axis, which is measured to 
be the zero field magnetic structire in CeCuBi 2 . J, J' and J±_ 
indicates three symmetry distinct spin exchange interactions 
between Ce magnetic moments. 


Magnetic susceptibility measurements on CeAgBi 2 
again show the presence of magnetic anisotropy but 
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smaller than in CeCuBi 2 , which indicates weaker Ising 
anisotropy of the localized moments. Furthermore, when 
Ag substitutes for Cu on the the transition metal site, 
spin superexchange between Ce moments occurs via the 
4d or 5s Ag orbitals, as compared to the 3d or 4s electrons 
of Cu, thereby changing the effective superexchange cou¬ 
plings J, J', Jj_. Moreover, the slight difference in lattice 
parameters and bond angles may also lead to changes in 
the RKKY interactions induced by coupling to the itin¬ 
erant electrons of Bi. In general, this confluence of com¬ 
peting effects may stabilize very complicated magnetic 
ordering. We also note that further-neighbor interactions 
besides J, J' and Jl generally induce incommensurate 
spiral ordering, affording an even richer set of possibili¬ 
ties. For the purposes of a preliminary analysis, we focus 
on our minimal spin model keeping just J, J' and J_l, and 
describing the possible ordered states as a function of the 
applied field. Notably, even this simple model exhibits a 
plethora of magnetic orders separated by a sequence of 
metamagnetic transitions. 

A convenient approach to understand the magnetic or¬ 
der is to separately consider every two layers of Ce sites, 
corresponding to a single unit cell lattice spacing along 
the c axis. These two layers of Ce sites then form a buck¬ 
led square lattice, where the nearest and next-nearest 
neighbor spin-spin couplings are J and J± respectively; 
we will refer to this as a ‘bilayer’. Reasoning in analogy 
with the well-studied Ji — J 2 square-lattice Heisenberg 
model, we anticipate that FM, Neel and stripy phases 
can be stabilized in each bilayer depending on the rela¬ 
tive signs and magnitudes of J and J_l . We may then use 
these orders as a building block to construct the full three 
dimensional magnetic order by coupling the Ce bilayers 
via an interlayer coupling J'. We now briefly summarize 
how competing interactions and magnetic anisotropy can 
lead to several metamagnetic transitions within this ap¬ 
proach, even if we restrict ourselves to working within 
a classical spin approximation. In order to have several 
metamagnetic phases, we require that J, J', J± and A 
are all of the same order; this is also consistent with 
expectations based on experiments and a study of lat¬ 
tice parameters. To describe the ordered states stabilized 
with this parameter regime, it is convenient to work with 
a magnetic unit cell 8 times larger than the original unit 
cell (i.e., containing 16 Ce sites), and optimize the clas¬ 
sical energy of this unit cell for a given set of parameters 
Jij, A and Ii, we minimize the energy of Eq.(l). 

Fig. 6 shows the magnetization m (m = 1 is the fully 
polarized classical spin per site) as a function of c-axis 
magnetic field for specific values of parameters J = — 1, 
J' = 2, J±_ = 1.5 and A = 1, which are the parameters 
that best describe the experimental data. In fact, we find 
six distinct phases over our field range. Note that the ex¬ 
act same argument can be applied when J and J' values 
are interchanged, due to site connectivity. The magnetic 
ordering pattern in each of these phases is illustrated in 
Fig. 7. 

At a zero field, each buckled square lattice bilayer 


forms a stripy phase, as | J| < | J_l| and the stripy phase 
is robust under a small and finite magnetic field. With 
the chosen parameters, this phase is degenerate in energy 
with the “up-up-down-down” like ordering (discussed for 
the case of CeCuBi 2 and illustrated in Fig. 5). Further 
increase of magnetic field, however, induces a transition 
into a phase where spins in the lower two layers are polar¬ 
ized along the field direction and spins in the upper two 
layers are anti-aligned to the field but have a finite cant¬ 
ing. The canted moments induce a finite magnetization, 
preventing perfect cancellation of magnetic moments be¬ 
tween the upper and lower bilayers. At « 4, a partial 
spin-flop transition occurs, yielding a finite magnetiza¬ 
tion jump. In this phase, spins in the upper two layers 
almost form a stripy phase, similar to the phase (I). In 
the lower two layers, however, all spins attempt to align 
along the field direction but four spins form a different 
canting angle with respect to the other four. By smoothly 
changing those canting angles, this phase reaches a half¬ 
magnetization plateau region, where spins in the upper 
bilayer form a stripy phase and spins in the lower bilayer 
are fully polarized. With stronger magnetic fields, there 
is another partial spin-flop transition exhibiting a finite 
jump in the magnetization. In phase (V), all 16 spins 
attempt to be aligned along the field direction but four 
of them have a different canting angle until their canting 
angle becomes zero and a fully polarized spin state (VI) 
is realized at a saturation field H « 14. 

Based on the above analysis, we conclude that com¬ 
peting exchange interactions and anisotropy can in prin¬ 
ciple lead to several metamagnetic transitions with field. 
We exemplified this for a specific parameter regime in 
Eq.(l) and found six distinct phases that may emerge 
taking into account the most dominant exchange interac- 



FIG. 6: Plot of magnetization m vs field h applied along c axis. 
In the presence of competing spin exchange interactions and 
anisotropy, there exist several magnetic phases with applying 
fields. (In particular, we take J = —1, J' = 2, J± = 1.5 and 
A = 1 in Eq.(l)) Fig. 7 describes different magnetic ordering 
patterns in every phase, I through VI. Inset shows experimen¬ 
tal phase diagram from Fig. 4. 















tions. Future neutron scattering studies will be valuable 
to confirm such scenario, and lay a firmer foundation for 
theoretical studies. More detailed theoretical and exper¬ 
imental analyses will be discussed in Ref. 16, which will 
also discuss scenarios beyond this minimal model and at¬ 
tempt to connect the magnetic ordering with transport. 

IV. CONCLUSIONS 

In conclusion, we have shown that high-quality 
CeAgBi 2 single crystals (Tjv = 6.4 K) present a rich field- 
tuned phase diagram with five metamagnetic transitions 
at 40 mK. In contrast to other CeTX 2 members, a 
strongly first-order transition is observed at ~ 79 kOe 
in the vicinity of the transition to the paramagnetic 
state. Remarkably, we unveil an anomalous Hall contri¬ 
bution for fields < 54 kOe which is drastically altered 


when H is tuned through a trio of transitions at 57, 78, 
and 84 kOe, suggesting that the Fermi surface is recon¬ 
structed in a subset of the metamagnetic transitions. Our 
results shed light on hidden properties of CeAgBi 2 and 
open new avenues for both experimental and theoretical 
studies on non-collinear magnetic structures, quantum 
(bi-)criticality, and Fermi surface effects, to name a few. 
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FIG. 7 : Magnetic ordering patterns for the phases (I) through (VI) in Fig. 6. For a clear description of ordering 
patterns, three figures are shown for every phase: the figure at left shows the ordering pattern in the 16-site magnetic unit cell, 
while the two figures on the right are top-down views of upper and lower Ce bilayers, respectively. (I) Spins in every bilayer 
form a stripy phase on a buckled square lattice. (II) Spins in the lower bilayer are polarized along the field direction while 
spins in the upper bilayer cant in the field direction. A small canting induces the finite magnetization m as seen in Fig. 7 . (Ill) 
Spins in the upper bilayer form a stripy phase and four spins flop in the lower bilayer, resulting in a partial spin-flop transition 
with a finite magnetization jump at the transition point h ~ 4. (IV, Half magnetization plateau region) Spins in the lower 
bilayer are polarized along the field direction and spins in the upper bilayer form a stripy phase. (V) Spins in the lower bilayer 
are polarized and four spins in the upper bilayer flop from a stripy phase, resulting in another partial spin-flop transition and 
magnetization jump near /i « 11. (VI) All spins are fully polarized along the field direction beyond the saturation field. 

































































